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1. Introduction 

In recent years there has been a resurgence of interest in charmonium spectroscopy, mainly 
due to the unexpected, narrow charmonium states discovered by BaBar, Belle and CLEO collabo- 
rations. The properties of these new states, called X, Y and Z, have been discussed in some recent 
reviews [1,2]. With the experiments at LHC and BES-III and the planned PANDA experiment at 
the new FAIR accelerator, more information about these states and the discovery of further states 
can be expected. Theoretical predictions, especially those from lattice QCD calculations, will pro- 
vide important guidance for these experimental efforts. 

Recent lattice studies of charm spectroscopy include [3, 4] and progress was reported by a 
number of groups at this conference [5, 6, 7]. In this work, we study the charmonium spectrum, 
including radially and orbitally excited conventional cc states and hybrids, on an anisotropic lattice 
ensemble generated by the Hadron Spectrum Collaboration [8]. We also study the contribution 
of disconnected diagrams. The distillation quark- field smearing algorithm [9], which is essential 
to improve the precision and enable efficient calculations of a broad range of hadron coiTclation 
functions, is used for the first time in charmonium spectrum study. By the application of the 
variational method [10] with a large basis of derivative-based interpolating operators, we aie able 
to extract and identify high spin (up to spin-4) states as well as exotic and non-exotic hybrids. 

This paper is organized as follows: we first give the details of the lattice setup and the tuning 
of parameters in the quark action, then explain the construction of the interpolating operators and 
the variational method as well as the method used to identify the extracted states. In section 2 
our results are presented by Z^*- while section 3 summarizes our calculation and describes ongoing 
work. 

1.1 Lattice Details 

To study high-spin states, exotic states and radial excitations the anisotropic lattice has proven 
very effective [11]. By making the lattice spacing in the temporal direction, at finer than the spatial, 
a.v additional resolution is gained in the direction in which measurements are made for resonable 
computational cost. In this work we use the lattices and parameters of the Hadron Spectrum Collab- 
oration which are described in detail in Ref. [8]. In brief, we use a Symanzik-improved gauge action 
and a tree-level tadpole-improved Sheikholeslami-Wohlert fermion action with stout-smeared [12] 
spatial links. The anisotropy t, = as/ui = 3.5 and the charm quark is simulated with the same 
(relativistic) action as used for the light quarks. Distillation [9] is useful in this study since it im- 
proves the signal-to-noise in traditionally noisy states such as hybrids and disconnected correlators 
and simplifies correlator constuction allowing for a large basis of operators with inversions on all 
timeslices. Some relevant parameters ai^e summarised in Table 1. 

Nf Volume A^cfgs m;r (MeV) t, A^ev ai\Q.) (G&N) 
^+1 16^ X 128 96 396 3^ 64 5.667 

Table 1: The reader is referred to Ref. [8] for more details. In this study we use 64 distillation eigenvectors 
and inversion on all timeslices (A'r = 128) in the distillation space. 
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1.2 Heavy quarks on an anisotropic lattice 

In dynamical simulations on anistropic lattices the bare anisotropy must be tuned, in the gauge 
and fermion sectors simultaneously, to ensure that the measured quantity in a physical probe takes 
its desired value; 3.5 in this case. This nonperturbative tuning by the Hadron Spectrum Collabo- 
ration is described in Ref. [13] where the dispersion relation for light quark masses was used to 
tune the anisotropy in the fermion action. We have investigated the accuracy of this bare parame- 
ter, tuned at mught for charmonium (and D) physics and find that the parameter tuning in the light 
sector leads to an underestimate of the slope in the charmonium dispersion relation of ~15%. Fig- 
ure 1 shows two different dispersion relations: for the rjc on the left and the D meson on the right. 
The left plot compares the dispersion relations of the rjc when the heavy quark perambulators are 
calculated with a bare anisotropy tuned at mught and at mcharm- 




Figure 1: The dispersion relations for the 77^- (left plot) and D (right plot) mesons. The left plot shows two 
dispersion relations: the points in red are determined using anisotropy parameters tuned at light quark masses 
while the blue points show the dispersion relation determined after a retuning of the fermion anisotropy at 
the charm quark mass. The measured anisotropics are also shown, determined from the slope by a linear fit 
including all points out to n^ = 4. The right plot shows the dispersion relation for the D meson using the 
parameters tuned at mcharm. on a subset of 15 configurations. The tuning of the anisotropy at mcharm also 
results in a small mass retuning which has been taken into account in these plots. 

K 

The measured values are shown on the plot, determined from Vf = — —, where b is the slope 

of the dispersion relation. The "retuned" bare anisotropy is then used to determine the dispersion 
relation for the D meson. We find the anisotropy measured from the slope is 3.38(3), in reasonable 
agreement with the target of 3.5. In passing we also note that the mass of the D meson is determined 
to be 1883(64) MeV (for unphysically-heavy light quarks). The measured anisotropy is precisely 
determined in both cases, with a statistical error less than a percent and excellent fits out to momenta 
of = 4. A more detailed analysis of the heavy-light meson spectrum is underway but some 
preliminary results on just 15 configurations, include the mass difference A{Mo* — M/))=137(4) 
MeV to be compared to the experimental result of 142.1 2(7 )MeV. 

1.3 Disconnected diagrams 

In this study we present our preliminary results for the disconnected contributions to the rjc, 
on a subset of 39 configurations. The rjc is a colour singlet and Wick contractions produce bubble 
diagrams, as well as the familiar connected diagrams, which should be included for an accurate 
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determination of the mass. It can be argued that in charmonium these disconnected diagrams are 
OZI suppressed, however to achieve a precision on mass determinations at the percent level it is 
valuable to determine their contribution more reliably. 

Distillation simplifies the calculation of disconnected diagrams considerably. The correlation 
function becomes a product of traces over rank M matrices, where M is the number of distillation 
eigenvectors used in the calculation. 



C^'''{t',t) = Tr [<I>^(t)T(t,t)] Tr [<I>^(t')T^(t',t')] , 



(1.1) 



where t(?,?) are the perambulators encoding the quark dynamics and the contain information 
about the spin structure of the meson. In Figure 2 we show the disconnected correlators which are 
relevant for the r\c as well as the connected coiTclator. Note that clear signals for the disconnected 
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Figure 2: The disconnected contributions to the r\c from the charm, strange and light sectors. The connected 
correlator is also shown. The disconnected correlators are determined in this preliminary work on just 39 
configurations but nevertheless a clear signal is determined which persists over a number of timeslices. We 
are currently increasing statistics for a full analysis. 

contributions from the charm, strange and light sectors all persist over a reasonable number of 
timeslices. One can write a matrix in flavour space to describe the mixing of these contributions 
and subsequently extract the full T]f correlator including all disconnected effects as well as the 
mixing angles, and this is ongoing work. 



1.4 Construction of interpolating operators 

The hadron spectral information is extracted from the time dependence of two-point coitc- 
lation functions between creation and annihilation operators. The simplest meson interpolating 
operators are the local fermion bilinears \lf{x)T'[if{x), where the quantum numbers are determined 
by the choice of gamma matrix, Y. However, these local operators allow for only a limited set 
of quantum numbers and to study higher spin and exotic states one must use non-local operators. 
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In this work, we use the derivative-based operators [3, 11, 14] which are constructed by using 
spatially-directed gauge-invariant derivatives within a fermion bilinear. The general form is 

Yt'^^'^j-Y, (1-2) 

^ y ^ y 

where D =D — D . The use of the "forward-backward" derivative is helpful to ensure that the 
operators have definite charge-conjugation at finite momentum, but not necessary for zero momen- 
tum which is the only case we consider in this work. It is straightforwai^d to build this type of 
operator with definite spin, parity and charge-conjugation in continuum space. The derivatives can 
be expressed in a basis transforming like spin-1: 

£>m=-l = ^(£>.v-?Z)y), D,„=o = lD„ D,n=l=--j={D , + lDy). (1.3) 

Using the standard SO (3) Clebsch-Gordan coefficients, one can combine any number of derivatives 
with gamma matrices to construct operators with the desired quantum numbers. The choice of 
gamma matrix plays a role in setting the parity and charge-conjugation and the relevant matrices 
for this study are listed in Table 2 with a naming scheme, as used in Ref. [3]. With one derivative. 
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Table 2: Gamma matrix naming scheme, 
the operators of spin {J,M) can be formed as 

(r X D^jliY-^ = X {Jr,mr; 1 , m2\J,M) Wmr'^rm Y, (1-4) 

where / is the continuum spin and M is the number of components of that spin. For scalar-like 
gamma matrices, 7r = 0,nir = 0. For vector-like gamma matrices, Jr = l,mr = —1,0, 1. Two 
derivatives can be combined as a definite spin Jj), and then combined with the gamma matrix as 

{TxD\iy^'^= {Jr,mr;JD,mD\J,M){\,mu\,m2\JD,MD)WmrD m2¥- (1-5) 

mY.mo,mi ,m2 

At the three-derivative level, we follow the convention in Ref. [11] for the ordering in which we 
couple the derivatives. The first and third derivatives are coupled to a definite spin /n and then the 
second derivative is coupled with /13 to form a definite spin Jp. 

{rxD^^lj^Y''' = £ {Jr,mr;JD,mD\J,M){l,m2;Ji3,mu\JD,MD) (1.6) 

^ y ^ y. ^ y 

X (l,mi;l,m3|7i3,Mi3)v/rmj.D„„ Dm^Dm^W- 

The operators built in this simple way have definite charge-conjugation. This is possible because 
the forward-backward derivatives are used. This procedure can be extended to any number of 
derivatives. In this work, we use operators with all combinations of up to three derivatives. 

In lattice calculations, states are not classified according to the spin (J, M), but rather according 
to the irreducible representations (irreps) of the cubic group. Different continuum spin components 
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Irreps(A) Ai A2 Ti T2 E 

dA I 13 3 2 

/ 0,4, 6, ••• 3, 6,7, ••• 1,3,4, ••• 2, 3,4, ••• 2,4,5, ••• 



Table 3: The table shows the irreducible representations A of the cubic group, together with their dimen- 
sions (c/a) and continuum spin content J. 

can be distributed across several irreps. Table 3 shows the five irreps, A of the cubic group and their 
continuum spin content. The continuum operators described above must then be subduced into 
lattice iiTcps. The subduction can be performed using linear combinations of the M components 
for each J: 

^AA-I^A>''''' (1-7) 
M 

where A = 1, • • • ,dA, is the "row" of the irrep A, ^■''^ is the continuum operator of spin {J,M) and 
are the subduction coefficients which form an orthogonal matrix Lm •^a^'^a'^I' ~ ^a,a'^a,A' ■ 
For the details of the derivation of the subduction coefficients see Ref. [11]. 

1.5 Variational method 

The variational method [10] is a very powerful tool to extract multiple excited states, utilizing 
a large basis of operators within a given symmetry channel. Mathematically this method involves 
solving the generalized eigenvalue problem 

C{t)Va = Xa{t)C{to)Va, (1.8) 

where C{t) is the correlator matrix at time slice t, Cij{t) = {iffi{t)S'j{0)). The eigenvectors ai^e 
orthogonal on the metric C{to), vlC{to)vh = 5a.b- The eigenvalues Xa{t) satisfy 

= (1 + ^(g-A'"<,(r-'o))^ (•1_9) 

where ni„ is the mass of a state labeled by a and Am^ is the distance of m„ from the mass of 
the nearest state. At large t, the principal correlators tend to ~ e^™a('^'o) practice, a second 
exponential is used in the fit form: 

Xa{t) = (1 -An)e-"'"^'-"'^ +A„e-'"«('-'°\ (1.10) 

where A„, rtia and m', are fitting parameters. 

Choosing a proper ?o plays an important role because we use a finite number of operators. The 
eigenvectors are forced to be orthogonal on the metric C(fo) by the solution procedure. This will 
only be a good approximation to the true orthogonality if the correlator at ?o is dominated by the 
lightest dim(C) states. In this work, we follow the "reconstruction" procedure described in Ref. [3] 
to choose ?o: the coiTclator matrix is reconstructed from the extracted masses and the eigenvectors, 
a ;^^-like quantity is defined and fo is chosen so that this is minimised. 
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1.6 Spin identification 

As the lattice spacing approaches zero, the rotational symmetry is restored and the spins can 
be identified by the emergence of energy degeneracies between different irreps in the continuum 
limit. For example, a spin-2 state would be expected to appear as degenerate energies within the T2 
and E irreps. However, this method requires repeating the calculation at a series of different lattice 
spacings. In this work we only have a single lattice spacing. 

An alternative approach valid for finite lattice spacing is to consider the "overlap" of the state 
extracted in the lattice irreps onto the operators, (0|(^|?i). From Eq. 1.7 one can see that the opera- 
tors constructed on the lattice carry a "memory" of the continuum spin from which it is subduced. 
If our lattice is reasonably close to the continuum limit, we can expect that an operator subduced 
from spin / only only overlap strongly onto those states of continuum spin J. Figure 3 shows the 
overlap histograms of four low-lying states in channel. It is clear that the first two states over- 




0.5370(0.0002) 0.6735(0.0011) 0.6774(0.0014) 0.7791(0.0039) 



Figure 3: The overlaps of four lowest-lying states onto the 26 different operators used in the Tj analysis. 
The color-coding indicates the continuum spin J from which the operators are subduced (red: J=l, blue; 
J=3, yellow: J=4). 

lap strongly onto 7 = 1 operators, the second state overlaps onto 7 = 3 operators and the fourth 
state overlaps onto 7 = 4 operators. 

To be more quantitative, we can compare the overlaps obtained in different irreps. The contin- 
uum operators have definite spin (0|^'''^|7',M') = Z^''^5jj'5m,m'- Therefore, as continuum sym- 
metries are restored lattice operators should obey (0|i^^'^|7',M') = s'j^^Z^^^Sj f. The value ofZ^ 
is common in different irreps. This suggests that we can identify a state by comparing the Z- values 
obtained in different irreps. For example, a 7 = 2 state created by 7 = 2 operator will have the same 
Z- value in each of the T2 and E irreps. 
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2. Spectrum Results 

The masses of charmonium states are calculated in every A^*- channel. Continuum spins are 
assigned using the method described above. In the following we will discuss each PC combination 
separately. The continuum spins are consistently colour-coded in the discussion that follows: black 
for 7 = 0, red for / = 1, green for J = 2, blue for 7 = 3 and yellow for 7 = 4. We use brown 
when we cannot make an unambiguous spin assignment. Figures 4 to 7 show our results in each 
PC combination, ordered by lattice iiTcp and colour-coded by continuum spin, J. The lowest 
thresholds for open-charm, DD and D^Ds are also shown. These thresholds have been calculated 
on the same lattices, described in Table 1. 



2.1 /++ 

Figure 4 shows the extracted 7++ states in lattice irreps. The lowest band of states, located 
4200 I 
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Figure 4: The plot shows the masses of the extracted 7++ states by lattice irreps. 



in the irreps Ai, Ti, Tj and E, are identified as the near-degenerate P-wave triplet 0^, 1^ and 
2++ by looking the overlaps of these states with the operators in each channel. The assignment of 
2++ is alsi 
channels: 

(ground state) 



2++ is also supported by comparing the overlaps onto the operator (p x T^yli)''"^ in the T2 and E 



Ztj (ground state) 



0.997(4). (2.1) 



This result clearly indicates that the ground states in T2 and E are different components of a 2++ 
state. The second band of states, at around 3950 MeV, are also located in the iiTcps Ai, Ti, T2 and 
E. At this energy range they could all belong to a single 4++ state considering the continuum spin 
contents of these four iiTcps. However, the overlap histograms indicate that they ai^e likely to be 
0++, 1++ and 2++ states. The 2++ assignment is also supported by comparing the overlaps onto 
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the operator (p x Oy^J-'"^ in T2 and E channel: 

Zf (1st exited state) , 

i (p X ^ , y-^) = 0.92(8). (2.2) 

Zr2(lst exited state) ^^'^ •'=^^ ; v y 

Another issue arising here is that this 2++ state could be either a radially exited P-wave or an F- 
wave. Checking the overlaps of this state onto an F-wave operator (p x ^/||=2 7=3)"'"^ -^2 and 
E, we find they are small in comparison to the overlaps onto the operator (p x D^J^j)''"^. Thus, 
identifying this 2++ state as a radially exited P-wave is preferable. Therefore, we can identify these 
near-degenerate 2++ states as the radially exited P-wave triplet. 

The remaining state in the A 1 channel can be identified as a 4++ state by its large overlap on 
the J = 4 operator (p x Dfj^^2 j=3y^^- The only state in A2 is probably a 3++. The remaining two 
states in Ti are expected to be a 3++ and a 4++. We can assign the second exited state as a 3++ and 
the third exited state as a 4++ by simply looking the overlap histogram in this channel. Similarly, 
the remaining three states in T2 channel are identified as a 2++, a 3++ and a 4++. The remaining 
two states in E are identified as a 2++ and a 4++. These assignments are all well supported by 
comparing the overlaps between different irreps. Comparing the overlaps of the assigned 2++ states 
in T2 and E onto the operator (p x £)|||^^2 7=3)^"^' S^*- i^stio to be 1.12(5). Comparing the 
overlaps of the assigned 3++ states in A2, Ti and T2 channels onto the operator (p x ^7^3=2 /=3)"'"^ 



, we have 



(ground state) nnAti\ '^^^ (ground state) 



0.94(2), ^ 7^^'" — — ^=0.99(2). (2.3) 



Zj-, (2nd exited state) ' Zj^ (3nd exited state 



For the 4++ state distributed in A\, T\, T2 and E channels, the overlaps onto the operator (p x 
^713=2 7=3)^"^ all appear to be compatible with a common state assignment: 

Za, (2nd exited state) Za, (2nd exited state) _ ^ Za^ (2nd exited state) _ ^ Q^/gN 

Ztj (3rd exited state) ~ ' Zj^ (4th exited state) ~ ' Z^ (2nd exited state) ~ ' 

(2.4) 

2.2 J 

Figure 5 shows the extracted / states organized by lattice irreps. The lowest state in the 
T\ channel can be identified as the J l\\f. The first excited state in T\ is also a 1 by looking the 
overlap histogram, which is very similar to the histogram for the ground state. Thus, we identify it 
as a radial exitation of the ground state, J /\if{lS). 

The band of states around 3850MeV, featuring the ground state in the A2 irrep which is safely 
identified as a 3 state, can be naturally explained as the D-wave triplet. The two very close 
states above 7/i//^(25') in the T\ channel ai^e expected to be a 1 and a 3 . The overlaps onto the 
operator (p x Df}^2y^^ gi^^ 

Za, (ground state) _ Za, (ground state) 



Z^i (2nd exited state) ' Z^, (3nd exited state) 

This shows that the third exited state should be a 3 while the second exited state is a 1 . The 
ground state in E, which is much lower than the other states in this channel, can be identified as a 
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Figure 5: The plot shows the masses of the extracted 7 states by lattice irreps. 

2 . The ground state in T2 is associated with the ground state in E by looking at the overlaps onto 

Zp (ground state) , , 

= 1-02 1 • (2-6) 

Zjj (ground state) 

The first exited state in T2 is associated with the ground state in A2 by the relation of their overlaps 

Za, (ground state) 



onto the 7 = 3 operator (p x D^^^j)"'"^- 



27-2(1 St exited state) 



1.04(1). (2.7) 



In the range above 4000 MeV, we identify three 1 states in the Ti channel by looking at the 
overlap histogram. The state with mass 4189(54) MeV, indicated in the plot with lighter shading, 
shows a dominant overlap onto the operator {n x £)yl])"'^^ which is built from a gluonic excitation 
coupled to a qq. Therefore, this state is probably not a conventional qq state, but could be a hybrid 
with non-exotic quantum numbers. As we will show below, possible non-exotic hybrids also appear 
in 0^+ and 2 ^ at the same energy range. Together with the exotic 1^+ state, they form a multiplet 
of near-degenerate hybrid mesons with J^^ = (0, 1,2) ^, 1 as discussed in Ref. [15]. The other 
two 1 states in Ti channel are possibly radial excitations of the J/Y- 

The overlap histogram of the Ai channel clearly indicates that the ground state is a 4 . We 
also find a 4 state in each of the Ti, T2, and E channels. Their overlaps onto a 7 = 4 operator 
{ai X £>lf3 ^2,7=3 gives the relations ZaJZt^ = 1.10(5), ZaJZt, = 0.93(3), ZaJZe = 1.03(3), 
showing that they are different components of a single 4 state. In this energy range, we would 
expect it to be a G-wave. 

2.3 J+- 

In the 7+^ sector, the 7-even states are exotic in the sense that they are not allowed in the quark 
model. In Figure 6 we show the extracted 7+^ states by lattice irreps. The lightest state in T\ is 
far lower than the other states and has no partners in other channels. We can safely assign it as the 
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Figure 6: The plot shows the masses of the extracted 7+ states by lattice irreps. 

1+^, the he. The set of states around 4000 MeV are assigned to be a 1+^, 3+^ pair. The ground 
state in A2 can be taken as a 3+^ considering the continuum spin content of A2. The ground state 
in T2 could be a 2+^ or a 3+^. Since 2+^ is exotic, we expect it be higher than a non-exotic 3+^ 
state. The overlap histogram also indicates that the ground state should be a 3+^. The first and 
second exited states in T\ channel are assigned as a 1+^ and a 3+^ respectively. These assignments 
are confirmed 
A2, Tl and T2. 



are confirmed by comparing the overlaps of the 3+ states onto a 7 = 3 operator b\ x D^j}^2y~^ 



Z,, (ground state) _ ^ ^^^^^^ Z, (ground state) ^ ^ ^^^^^ 



Zjj (2nd exited state) Zt-, (ground state) 

We assign the ground state in Ai to be an exotic 0+^ state. The first exited state in T2 are associated 
with the ground state in E by overlap analysis. They complete an exotic 2+^ state. 

2.4 J-+ 

In the 7~+ sector, the 7-odd states are exotic. Figure 7 shows the extracted 7^+ states by lattice 
irreps. The lowest two states below the DD threshold in the A\ channel can be identified as 0^+, 
they are both S-wave, one is T]f (15) and the other one is 77^(25). The second exited state is also a 
^ , probably another radially excited S-wave according to the similarity of its overlap histogram 
with the lowest two states. The third exited state has large overlap onto the operator (p x D^jL^y^^ 
which contains gluonic excitation so this could be a hybrid state with non-exotic quantum numbers. 

The ground state in T2 and E can be easily identified as a 2^+ D-wave. The first exited 
states in these two channels are also assigned as 2 ^ and in this energy range could be a radially 
excited D-wave 2^+ state. However, looking carefully at the overlaps, we find that they both have 
big overlaps onto the operator (p x D^J^^^Y^^. Furthermore, the ratio of their overlaps onto this 
operator is Zj^IZe = 1.2(1) suggesting that it is probably a 2^+ hybrid state. 

The ground state in Ti is assigned to be an exotic 1^+ state. A non-exotic 4 ^ state is located 
very closely. In Ai, T2 and E channels, we also find a 4^+ state. Their overlaps onto the 7 = 4 
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Figure 7: The plot shows the masses of the extracted 7 + states by lattice irreps. 

operator {bi x ^7^3=2 7=3)''"'* §i^^ relations: 

1^ = 1.07(4), |ii = 1.09(5), ^=0.96(3), (2.9) 

Lt^ Le 

indicating that they are different components of a 4^+ state. The discrepancy of the mass in differ- 
ent irreps is possibly a discretisation effect. 

The ground state in A2 could be an exotic 3^+ state or a non-exotic 6^+ state. Since we do not 
use operators with overlap onto 7 = 6 states, a 3^+ interpretation is preferable. We also find a 3^+ 
state in T\ channel, which is associated with the ground state in A2 by overlap analysis. A 3^+ is 
also expected in T2 channel to complete this state. However, we don't have convincing evidence to 
unambiguously identify a 3^+ in T2. 



3. Summary 

We have presented a progress report on our study of the charmonium spectrum using anisotropic 
lattices and 2 + 1 fermion flavours. The spin-identification, operator construction and variational 
analysis used here is crucial for an accurate identification of states. In particular we are able to 
reliably identify spin 4 states and both exotic and non-exotic hybrids. 

We have investigated the dependence of the bare fermion anistropy on the quark mass and 
find that a re-tuning for heavy quarks is necessary. The anisotropy measured from the charmonium 
dispersion relation is then 3.50(2) for a target anisotropy of 3.5. It is also interesting to note that 
the dispersion relation is fully consistent with a relativistic interpretation and momenta up to n = 2 
can easily be included in the fit. 

Figure 8 is a summary of the spectrum discussed in the preceeding sections. We note that a 
number of hybrid states, exotic and non-exotic, are determined and we find strong evidence for a 
charmonium "super-multiplet" as discussed in Ref [15] for light mesons. The statistical resolution 
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Figure 8: Summary of the all the charmonium and exotic states. The dashed purple lines indicates the DD 
and DsDs thresholds. The red bars are the experimental values. 

of our determination of these states, with 96 configurations, is ~ 17 MeV. For states below threshold 
the precision is of the order 1% or less: ~ 1 MeV on the r]^ and J /\lf. 

A complete description of all disconnected effects in the Tj^ is also outlined above. Using 
distillation this calculation is considerably simplified and signals persisting over 5-10 timeslices 
on just 39 configurations are resolved. This ongoing work with higher statistics will be more 
completely described in a further paper. 
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